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Abstract. We introduce a general modeling framework to predict the outcomes, 
at the population level, of individual psychology and behavior. The framework 
prescribes that researchers build a cost function that embodies knowledge of what 
trait values (opinions, behaviors, etc.) are favored by individual interactions under 
given social conditions. Predictions at the population level are then drawn using 
methods from statistical mechanics, a branch of theoretical physics born to link 
the microscopic and macroscopic behavior of physical systems. We demonstrate our 
approach building a model of cultural contact between two cultures (e.g., immigra- 
tion), showing that it is possible to make predictions about how contact changes the 
two cultures. 
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1. Introduction 



Modern society features an increasing degree of interaction between 
cultures ("cultural contact") owing to, e.g., communication technolo- 
gies, immigration and other socio-political forces (Lull, 2000). In many 
countries cultural contact is perceived as both an opportunity and a 
threat and is related to important public issues such as immigration 
management and the need to "protect" national culture (Cornelius 
et al., 1994; Givens and Luedtke, 2005). Our understanding of these 
phenomena is, however, limited: we cannot predict the outcome of 
cultural contact, nor make plausible conjectures that can be used in 
policy making. Within this context, the aim of this paper is twofold: 
we first describe a general mathematical framework for modeling social 
interactions, then we make specific assumptions relevant to studying 
immigration, i.e., social contact between two groups that, typically, 
differ both in culture and relative size. 
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2. General framework 



2.1. A FORMALIZATION OF SOCIAL INTERACTIONS 

For simplicity, we focus on a single cultural trait, which may represent 
an idea, an opinion or a behavior that has two mutually exclusive forms. 
A useful example to keep in mind is being in favor or against an issue 
such as the death penalty, or any other issue that might be the subject 
of a yes/no referendum vote. Our framework allows to consider multiple 
traits without conceptual differences, although model analysis may in 
general be much more difficult. 

We consider a population of N individuals, labeled by an index i 
ranging from 1 to N. We associate to the i-th individual the variable 
Si, which may take values +1 or —1 representing the two possible trait 
values. For instance, +1 might represent a yes vote in a referendum, 
and —1 a no vote. The state of the whole population is thus encoded 
in an array s of N numbers, such as s = {+1, — 1, +1, . . .}. 

The hallmark of social interactions is that individuals may change 
their opinions or behavior owing to interactions with others. A given 
couple (i, k) can in principle be in one of the four states {+1, +1}, 
{— 1,+1}, {+1,-1} and {—1,-1}, but these outcomes, in general, do 
not have the same probability. Which one is more likely will depend on 
the characteristics of individuals such as their culture and personality. 
Our starting assumption is that individuals have no particular bias 
towards +1 or —1 opinions: what matters most is whether, by adopting 
one or the other value, an individual is in agreement or disagreement 
with others. There are two reasons for this assumptions. First, social 
psychologists have shown that, in most cultures, agreement or dis- 
agreement with others is a powerful determinant of individual opinions 
and behavior, often more important than holding a particular opinion 
(Bond and Smith, 1996); we will expand on this point in our model of 
immigration below. Second, our framework allows to introduce biases 
that favor a particular trait value, if needed. Indeed, any model in 
which individuals are biased can be recast as a model with unbiased 
individuals, plus an additional "force" that orients individual opinions. 
Thus our starting assumption of unbiased individuals does not reduce 
the generality of the framework. Again, we will make a specific example 
for the case of immigration below. 

To formalize these notions, we assume that individuals take on 
the trait that minimizes a cost function. We define the cost for 
individual i to agree or disagree with individual k as 
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where Jik is a number that summarizes the nature of the interaction 
between i and k, as follows. When i and k agree (siSk = 1) we have 
a cost = —Jik, while when i and k disagree (sjSfc = —1) we have 
Hik = Jik- Thus whether agreement or disagreement carries the lesser 
cost depends on the sign of Jik : Jik > favors agreement while Jik < 
favors disagreement. The magnitude of Jik gives how important it is for 
i to agree or disagree with j. If, for instance, Jik > Jim then it is more 
important for i to agree with k than with m, while Jik = means that 
agreement with k is not relevant to i. The signs and magnitudes of the 
Jjfc's become important when we consider that an individual interacts 
with many others. In this case, we assume that the costs Hik relative 
to each interaction sum up to a total cost Hi 

Hi(si) = ^2H(si,s k ) = -^JikSiSk (2) 

k k 

As anticipated above, we can take into account additional factors that 
may influence individuals modifying equation (2) as follows: 

Hi(si) = -^2 JikSiSk - hsi (3) 
k 

meaning that individual % is subject to an additional "force" that fa- 
vors Si = 1 if hi > and Sj = — 1 if hi < 0. The quantity hi may 
represent any factor that is not explicitly taken into account by the di- 
rect interaction with other individuals. For instance, it may summarize 
the influence of media, government campaigns or existing culture (see 
below) . 

We can now write a population-level cost function as the sum of 
individual cost functions: 

H(s) = J2 H i( s i) = ~J2 JikSiSk ~J2 hiSi ( 4 ) 

i i,k i 

We stress that the cost function is a theoretical computational aid 
to track which trait values are favored by the interactions Jik and 
the external forces hi. We do not assume that individuals explicitly 
compute or are aware of such costs. Rather, H(s) should be designed 
so that its minima correspond to those trait values that are favored by 
the mechanisms with which individuals interact. 

2.2. The role of statistical mechanics 

Once a cost function has been specified, it is possible to calculate 
population level quantities such as the average trait value using the 
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methods of statistical mechanics, a branch of theoretical physics. Sta- 
tistical mechanics was originally concerned with deriving the laws of 
thermodynamics from the behaviour of atoms and molecules (Thomp- 
son, 1979), but can actually be applied to understand the macroscopic 
(population level) properties of any system composed of many parts 
that interact according to given microscopic (individual level) rules. 
More recently its methods have found application in fields as diverse as 
biology (Mezard et al., 1987), neuroscience (Amit, 1989; Arbib, 2003), 
economy and finance (Bouchaud and Potters, 2000) and also social 
science (Durlauf, 1999). The starting point is to assign to each system 
configuration s a probability Pr(s) according to the Boltzmann-Gibbs 
formula (Thompson, 1979) 

e -H(s) 

where the sum runs over all possible configurations of the system. 1 

By means of (5) a given configuration is considered more or less 
likely according to whether it is more or less costly: a low value of H(s) 
results in a high probability of s, and vice- versa. Assigning probabilities 
based on a given cost function is the heart of statistical mechanics and 
is inspired by the principles of thermodynamics (see the appendix for 
a short discussion, and Thompson, 1979, for a fuller treatment). 

Once a probability is assigned to system configurations, it is possible 
to compute the expected values of quantities of interest and to relate 
them to the parameters that describe the system. For instance the 
average cultural trait defined by 

m ( S ) = JyH S i ( 6 ) 

i 

would have an expected value given by 

to = ^ m(s) Pr(s) (7) 

s 

Note that, while m(s) is the average trait value in a given configuration, 
to is the average trait value over all possible system configurations, 
each one weighted according to its probability of occurrence. These 
probabilities depend on the cost function H and thus on the parameters 
that appear in its expression, i.e., the Jjfc's and h^s. 

Rather than directly attempting to compute expected values such 
as (7), statistical mechanics aims to compute the so-called free energy 
of a system, defined as 

/ = -log (8) 
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The rationale for this strategy is that important quantities such as (7) 
can be easily computed from knowledge of the free energy function, typ- 
ically by taking suitable derivatives with respect to system parameters 
(see appendix). The basic task of statistical mechanics is thus, after the 
cost function H has been specified, to calculate the free energy function 
for a given system. We stress that the form of the cost function is not 
given by statistical mechanics; rather, it is the outcome of a modeling 
effort relative to a specific problem. We now make an example of how 
one may proceed. 



3. The cultural outcomes of immigration 

We illustrate here the potentials of our framework considering the 
impact of immigration on culture. We consider a large and a small 
population, which will be referred to, respectively, as residents (R) 
and immigrants (I). We let Ni be the number of residents, and N2 of 
immigrants, with Ni 3> N2 and N = N\ + N2 the total number of indi- 
viduals. We are interested in how cultural contact changes the average 
trait values in the two populations, with the aim of understanding the 
effect of one culture upon the other. 

3.1. Interactions between residents and immigrants 

Our main assumption regarding how residents and immigrants interact 
is that people, generally speaking, tend to agree with those who are 
perceived as similar to oneself and to disagree with those perceived as 
different. In social psychology this is known as the similarity-attraction 
hypothesis. It has received ample support, although the details of how 
we interact with others often depend on social context (Michinov and 
Monteil, 2002; Byrne, 1997). We consider this assumption a general 
guideline, and in modeling a specific case it can be modified without 
difficulty. We formalize the similarity-attraction hypothesis by assum- 
ing that high perceived similarity corresponds to positive values of 
Jik, and low perceived similarity to negative values. Since residents 
and immigrants have generally different cultures, we may assume the 
following structure for the interaction coefficients J^. We let the inter- 
action between any two residents be jC 1 ' 1 ) > 0; the similarity-attraction 
hypothesis suggests that this be a positive number, whose magnitude 
reflects how strongly residents prefer to agree among themselves. Like- 
wise, we let the interactions between immigrants be j( 2,2 ) > 0. The 
mutual interactions j( 1,2 ) and j( 2,1 ) should model whether residents 
prefer to agree or disagree with immigrants, and vice- versa, and how 
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strongly so. If resident and immigrant cultures are very different, the 
similarity-attraction hypothesis suggests to take both J^ 1 ' 2 ) and J^ 2 ' 1 ) 
as negative, but the best choice of values depends on the specific case 
one intends to model. 

Note that we are assuming that depends only on population 
membership and not on the particular individuals i and k (the so- 
called mean field assumption in statistical mechanics) . This assumption 
greatly simplifies mathematical analysis but is not wholly realistic. It 
can capture the average variation in interactions across population but 
not the variation that exists within each population. For instance, 
a more realistic assumption would be to take the JikS as random 
variables whose mean and variance depend on population member- 
ship. We plan to return on that model (which would represent the 
two-population generalization of the Sherrington-Kirkpatrick model in 
statistical mechanics, Mezard et al., 1987) in future studies. 

When modeling interactions, a technical requirement is that the 
value of the cost function be proportional to total population size N. 
This guarantees that the free-energy function and important quantities 
such as average trait value, equation (7), scale appropriately with N. 
In our case the appropriate scaling is 1/2N, hence the interactions are: 



Jik = ' 



2N 

J^ 2 ) 
2N 

j^i) 

2N 

J( 2 ' 2 ) 
27V 



i,k G R 
i G R,k G I 
i g I,k G R 
i,k G I 



(9) 



3.2. Modeling pre-existing culture 

Before the two populations start to interact, residents and immigrants 
are each characterized by a given average trait value, say rh\ and rfi2, 
respectively. We consider rh\ and rfi2 as experimental data about the 
beliefs or behavior of each population, which could be obtained from, 
say, a referendum vote on a particular issue (e.g., the death penalty) 
or from statistical sampling of the population. 

That a population is characterized by a given average value fhi ^ 
means that the two forms of the trait are not equally common. 
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Specifically, the individuals with the +1 form are N(l + mj)/2, while 
N(l — rh\)/2 individuals have the —1 form. Pre-existing culture, in 
other words, is like a bias or force that favors one trait value over the 
other. For modeling purposes, it is convenient to describe pre-existing 
culture as a "cultural force" that acts to orient the opinion of otherwise 
unbiased individuals. This is possible including a force term in the cost 
function, as shown in (3). By standard methods of statistical mechanics 
(see appendix) it is possible to show that the force term corresponding 
to a particular average opinion rhi is 

h® = atanh(mi) - J^m (10) 

where atanh is the inverse hyperbolic tangent function. To summarize, 
a model in which individuals are biased so that the average opinion is 
rhi is equivalent to a model with unbiased individuals subject to a force 
given by (10). 

3.3. Model analysis 

So far we have specified interaction terms to model cultural contact 
between two populations and we have introduced equation (10) to rep- 
resent the pre-existing culture in the two populations. The next step is 
to compute the average trait values m\ and m<i in the two populations 
after immigration has taken place. The same method that allows to 
derive equation (10) enables us to derive the following equations for 
mi and 771,2 (see appendix): 

mi = tanh ((1 - a)J^ 1) m 1 + a(J (1 ' 2) + J (2 ^)m 2 + (11a) 

1712 = tanh ((1 - a)(J (2 ' 1} + J^mi + aJ^m 2 + h&) (lib) 

where a = N 2 /N is the fraction of immigrants in the total population 
and tanh is the hyperbolic tangent function. The values of m\ and m 2 
predicted by (11) depend of course on values of the J and h parameters, 
and on a. We give here a qualitative description of the different regimes 
that one can observe varying these parameters. We refer to Cohen 
(1974) for a proof of the following statements, in the context of an 
analogous model from condensed matter physics. 

The two key parameters are a, the fraction of immigrants, and the 
overall scale of the interactions J^, which we label J. If J is below a 
critical value J*, equation (11) has always one pair of solutions, for all 
values of a. In this case the two populations are essentially merging 
into a homogeneous one, with average cultural trait in between the 
two initial ones — more toward one or the other according to the value 
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Figure 1. Possible outcomes of cultural contact. Average trait value in the two 
populations, m = (1 — a)mi + am2, is plotted as a function of a (fraction of 
immigrants) and J (scale of the interaction) for an interaction matrix of the form (9) 
with J ( - I ' fc ) = J. Before the interaction the two populations have initial trait values 
of mi = 0.5 and rfi2 = —0.5. Only the most likely outcome is plotted (see text). 

of a. This regime is not surprising and corresponds to the naive pre- 
diction that one could have made a priori without applying statistical 
mechanics. 

If the interaction scale is large (J > J*), however, model predictions 
are highly non-trivial, suggesting that the outcome of cultural contact 
can be surprising. Depending on J there are two critical values for 
a: a\{J) and a<2{J) that delimit qualitatively different behavior. For 
a < a\{J) the resident culture dominates dominant and the immigrant 
culture disappears, i.e., 7712 is close to mi irrespective of the initial value 
m<i- The converse happens when a > ct2(J), i.e., immigrant culture 
dominates. The most interesting case occurs when ai(J) < a < a2(J)- 
In this regime (11) has two distinct solutions in which either of the two 
cultures dominates. That is, both cultures may survive the immigration 
process, generally with a different probability determined by system 
parameters. 

The parameter values that favor the resident or immigrant culture, 
have still to be worked out and will be the topic of future work. 2 Here 
we analyze the case in which the intensity of the interactions is uniform 
both within and between groups, = J. This is interpreted as two 

groups that do not really discriminate between themselves, so that 
disagreement with any particular individual carries the same cost inde- 
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pendent of which group the individual belongs to. We assume, however, 
that the two groups have initially a very different average trait value: 
rh\ = 0.5 and rfi2 = —0.5. In figure 1 we explore this system by plotting 
the average trait value after the interaction, m = (l-a)mi+am2, for a 
between and 10% and for J > 0. 3 For J = (no interaction) m is sim- 
ply the weighted sum of pre-existing trait values, m = (1 — a)fhi + am2, 
where each group contributes according to its size. As a function of 
a, this is a straight line. As the interaction increases the line slowly 
bends and for higher values of a we see a slight exaggeration of the 
pre-existing opinion rh\ (the surface in figure 1 rises slightly over the 
level m = 0.5). When J crosses a critical value J* ~ 1.125, however, 
a dramatic phenomenon occurs: the population undergoes a sudden 
change of opinion, switching to a value of m that is closer to, and indeed 
exaggerates, the initial value in the immigrant population, rfi2 = —0.5. 
Note that this sudden transition occurs for all values of a, i.e., irre- 
spective of the proportion of immigrants. The solution with m closer 
to m\ is still available to the system (not plotted in Figure 1), but as 
J grows past J* it is less and less likely that the system remains in 
such state (technically, for J < J* the solution with m ~ rh\ has a 
higher free-energy than the solution with m ~ rfi2 and thus becomes 
metastable, allowing fluctuations to cause a transition between the two 
solutions). Thus, according to this model, to prevent dramatic changes 
in resident culture, it would do little to restrict immigration (the effect 
of a is small in the graph). Rather, one should concentrate in reducing 
the scale of the interaction J, i.e., the strength of attitudes within and 
between groups. 



4. Discussion 



Attempts to apply mathematical-physics methods to social sciences 
have appeared in the litterature since the pioneering work of Galam 
et al. (1982). In this paper we have focused on statistical mechanics as a 
tool to bridge the gap from individual-level psychology and behavior to 
population-level outcomes. Our framework prescribes that researchers 
build a cost function that embodies knowledge of what trait values 
(opinions, behaviors, etc.) are favored by individual interactions under 
given social conditions. The cost function, equation (4), is defined by 
a choice for the interactions Jjfc and the fields hi that represent so- 
cial forces influencing individual opinions and behavior. This modeling 
effort, of course, requires specific knowledge of the social issue to be 
modeled. After a cost function has been specified, the machinery of sta- 



QlQt070606.tex; 2/02/2008; 7:37; p. 9 



10 



P. Contucci, S. Ghirlanda 



tistical mechanics can be used to compute population-level quantities 
and study how they depend on system parameters. 

We have demonstrated our framework with an attempt to under- 
stand the possible outcomes of contact between two cultures. Even the 
simple case we studied in some detail the model suggests that cultural 
contact may have dramatic outcomes (figure 1). How to tailor our 
framework to specific cases, and what scenarios such models predict, is 
open to future research. 
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Notes 

1 In physical applications, the parameter (3 (inverse temperature) usually mul- 
tiplies H in (5). Here, as in other applications of statistical mechanics such as 
combinatorial optimization (Mezard et al., 1987), we include f3 in H itself as an 
overall scale factor. 

2 As noted above, Cohen (1974) has studied a formally identical model arising 
from a condensed-matter problem, but there the hi's were free parameters, while 
here they are determined in terms of the and the m; through (10). 

3 The maximum admissible value for J is given by the condition that (10) has a 
stable solution, corresponding to the assumption that each culture is in equilibrium 
before cultural contact. This condition yields J ;$ 1.3. 
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Appendix 

A. Model solution 

It is a standard result of statistical mechanics (Thompson, 1979) that 
the free energy function of a system defined by a cost function of the 
form 

ff (*) = -2# E^-^E* ( 12 ) 

i,k i 

is obtained for the value of m that minimizes the function 

. Jo, 1 + m, 1 + m 1 — m , 1 — m , , 

i 7( m) = __ m 2_ /lm + __ log __ + __ log __ (13) 

The minimization of this function with respect to m yields the con- 
dition (10) which relates m and h and the Hamiltonian parameters. 
The structure of the free energy (13) admits the standard statistical 
mechanics interpretation as a sum of two contributions: the internal 
energy (the average of the cost function) 

U(m) =77 = -^m 2 - hm (14) 

minus the entropy 

or \ \^t> i m t> t \ l + m 1+m I - m 1 - m 
5(m) = - 2^ Pr (s) log Pr(a) = — log — — log 



2 
(15) 
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One can indeed show that the distribution function (5) may be deduced 
from the second principle of thermodynamics i.e. as the distribution 
for which the entropy is minimum at an assigned value of the cost 
function (Thompson, 1979). Equation (11) is obtained similarly from 
the representation of the free energy of the two population system as 
the minimum of the function 



F(mi,m 2 ) = -(1 - a) 2 J^4 ~ a 2 J^4 



-a(l-a)(J^ + J( 2 ' 1 ))m 1 m 2 
— (1 — a)h\m\ — ah2rri2 
+(1 - a) [+ 1± §^ log + ^ log 
+ al+l±mz log iif^ + log ±=p] 



1— mi 
2 



(16) 



The minimum condition yields (11). 



Qiqt070606.tex; 2/02/2008; 7:37; p. 12 



